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Abstract — We have investigated the main scattering mecha- 
nisms affecting mobihty in graphene nanoribbons using detailed 
atomistic simulations. We have considered carrier scattering due 
to acoustic and optical phonons, edge roughness, single defects, 
and ionized impurities, and we have defined a methodology 
based on simulations of statistically meaningful ensembles of 
nanoribbon segments. Edge disorder heavily affects mobility at 
room temperature in narrower nanoribbons, whereas charged 
impurities and phonons are hardly the limiting factors. Results 
are favorably compared to the few experiments available in the 
literature. 

Keywords - Low-field mobility, graphene nanoribbons, 
scattering, edge roughness, defects, impurities, phonons. 

I. Introduction 

Two-dimensional (2D) graphene sheets have demonstrated 
really attractive electrical properties like high carrier mobil- 
ity HI, lEl and large coherence length [3|. However, experi- 
mental data of mobility available in the literature show huge 
dispersion, ranging from lO'^ to 10'' cm^/Vs at room tempera- 
ture, signaling that the fabrication process is still poorly opti- 
mized and not fully repeatable. To guide process optimization, 
an exhaustive interpretation of physical mechanisms limiting 
mobility would be extremely useful. For Graphene NanoRib- 
bons (GNRs) a comprehensive experimental characterization 
of mobility is still lacking, mainly due to the difficulty in 
patterning in a repeatable way very narrow ribbons. Few recent 
interesting experiments are reported in [4 | and [5 j. GNRs may 
also suffer significant degradation of mobility due to additional 
scattering mechanims, such as edge roughness. 

The single most important aspect that makes graphene 
interesting for nanoscale electronics is its very high mobility. 
It is therefore of paramount importance to understand if 
also nano structured graphene can preserve the high mobility 
(often) measured in graphene sheets, much larger than that 
of conventional semiconductors. In addition, one would need 
to understand the effect on mobility of different options for 
graphene functionalization, which could be required to open 
a semiconducting gap in graphene. 

In the current situation, theoretical investigations ||6], Q 
and numerical simulations |8|, |9|, |10| can represent a 
useful tool to assess the relative impact of different sources 
of non-idealities on mobility and consequently on device 
performance, to provide guidelines for the fabrication process 
and a realistic evaluation of the perspectives of graphene in 
nanoelectronics. 



An analytical method and a Monte Carlo approach have for 
example been adopted in order to study line-edge roughness 
(LER) and phonon scattering-limited mobility in Ref. ||6l and 
Ref. fl], respectively. However, due to the reduced width of the 
considered devices, effects at the atomistic scale are relevant, 
therefore accurate simulation approaches like semi-empirical 
tight-binding are needed. 

In this work we present atomistic simulations of GNR- 
FETs, considering GNR widths ranging from 1 to 10 nm, 
and including scattering due to LER, single defects, ionized 
impurities, acoustic and optical phonons. A direct comparison 
with recently fabricated devices 14] will also be performed. 
Statistical simulations performed on a large ensemble of 
nanoribbons with different occurrences of the spatial distri- 
bution of non-idealities show that phonons, LER and defects 
scattering can likely explain the few available experimental 
data |4|, where mobility is down to the level of mundane 
semiconductors (order of 10^-10'^ cm^/Vs). 

II. Methodology 

A long GNR-FET channel, where mobility is properly de- 
fined, is given by a series of N GNR segments of length L like 
those we have considered in the simulation (Fig. [TJ. For the 
i-th GNR segment, the resistance Ri — Vds /h is the sum of 
two contributions, the channel resistance Rch,i and the contact 
(ballistic) resistance Rb — Vds/Ib (Ri = Rch,i+RB), where 
Vds is the drain-to-source voltage, whereas li and Ib are 
the total current and the ballistic current in the i-th segment, 
respectively. Assuming phase coherence is lost at the interface 
between segments, the resistance Rtot of the long channel 
GNR is therefore the sum of N channel resistances and one 
contact resistance, i.e.: 



N 



Rt. 



= J2 ^-'^-^ ]+Rb = N{R) -(N-l) Rb 



(1) 



where (R) — {l/N) X^ili Ri is the mean resistance evaluated 
on the ensemble of nanoribbon segments. Therefore, the 
mobility of a long channel would read: 



fJ-n 
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(2) 



Qtot Qtot N{R) - (TV - 1) i?B ' 

where the index n denotes each type of scattering mecha- 
nism limiting mobility (defects, edge -roughness or impurities). 
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Fig. 1 . GNR-FET channel with length NL and the simulated GNR segment 
with length L. 



Ltot = NL is the total GNR length, Qtot = E^=i Qi = 
N (Q) is the total charge in the channel and (Q) is the mean 
mobile charge in a segment. 

For large values of N, one can discard 1 with respect to 
N in Eq. (|2]) so that we obtain the formula we use in the 
paper ifTTI : 

{{R)-Rb){Q)' 

The root mean square error of mobility has been computed 
by means of a Taylor expansion up to the first order of Eq. (|3) 
with respect to statistical fluctuations of the resistance R = 

Rch + Rb- 



dR 

and therefore 
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where Ai? = and 



N 
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l/(iV-l)^(i?,-(i?)f 



(4) 



(5) 



(6) 



is the variance of R. 

Statistical simulations of resistance on a large ensemble of 
nanoribbon segments with different actual distribution of non- 
idealities have been performed. In particular, the mobility /i„ 
has been computed in the linear transport regime, for large gate 
voltages (Vgs) and small drain-to-source bias Yds — 10 mV. 
Mobility has been extracted by means of Eq. (|3]l considering 
an ensemble of = 600 nanoribbon segments with different 
disorder realizations for 1.12 nm-wide GNRs. Due to the 
computational cost, at least 40 nanoribbons segments have 
been instead simulated for 10.10 nm-wide GNRs. 

Statistical simulations of random actual distributions of de- 
fects, LER and ionized impurities have been computed through 
the self-consistent solution of 3D Poisson and Schrodinger 
equations within the NEGF formalism, with a -p^ tight- 
binding Hamiltonian ifTOl , extensively exploiting our open- 
source simulator NanoTCAD ViDES |12|. In particular, we 
have imposed at both ends of the segments null Neumann 
boundary conditions on the potential, and open boundary 
conditions for the transport equation. 



In order to compute the LER-limited mobility [Iler, sta- 
tistical simulations have been performed considering a given 
fraction H of single vacancy defects at the edges, ii is defined 
as the probability for each carbon atom at the edges to be 
vacant. In practice, each sample of nanoribbon with edge 
disorder is randomly generated assuming that each carbon site 
at the edges has a probability H to be replaced by a vacancy. 
Null hopping parameter has been imposed in correspondence 
of a defect at the edge. 

Defects have been modeled using the on-site energy and 
the hopping parameter extracted from DFT calculations [13 1. 
In particular, for a fixed defect concentration rid, each sample 
of defected nanoribbon with defects is randomly generated 
assuming that each carbon atom has a probability n^, to be 
replaced by a vacancy. 

As previously assumed in ab-initio calculations 1141 . we 
have considered a surface impurity distribution of positive 
charges equal to -1-0.4 q placed at a distance of 0.2 nm 
from the GNR surface, where q is the elementary charge. 
Again, if jiimp is the impurity fraction, a sample with surface 
impurities is randomly generated by assuming that each carbon 
atom has a probability njM p to be at 0.2 nm from an impurity 
in the dielectric layer 

In Figs. |2^-b, we show the distributions of Q when consid- 
ering line-edge roughness {H = 5%) and defects {rid = 2.5%) 
for W — 1.12 nm. In each picture we show the mean value 
(Q) and the standard deviation <tq of the random variable 
Q. For comparison, the corresponding normal distribution is 
shown. 

Phonon-limited mobility fiph (both acoustic and optical) 
has been computed by means of a semi-analytical model as 
in |6|, but extending the Kubo-Greenwood formalism beyond 
the effective mass approximation and accounting for energy 
relaxation at GNR edges fTSl . Starting from the Boltzmann 
transport equation, the phonon limited mobility for a ID 
conductor can be espressed as [16] : 



1 



dkx f (kx) 



(7) 



where v^j = dEj/dk^ is the electron velocity in the 

longitudinal direction x for the j-th electron subband and Tpj 
is the corresponding momentum relaxation time for electron- 
phonon scattering. In Eq. {■) denotes the expectation value 
averaged on the Fermi factor / as: 



(.9) = 



+ 00 



dkx -^g{kx) f{kx) 



(8) 



where riin is the one dimensional (ID) carrier density. In 
order to compute Eq. O, the following electron dispersion 
curve has been exploited for the j-th subband ifTSl : 



Ej (kx) = JeI^„ 
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(9) 



where Ecj — Ecjo — Q^c is the cut-off energy of the j- 
th subband when the electrostatic channel potential is 
different from zero {Ecj — Ecjo for — V). According 
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w=1.12 nm 

<Q>/q= 0.8809 
OQ/q=0.0642 




Fig. 2. The distributions of charge Q in each GNR segment {W = 1.12 nm) 
obtained from statistical simulations of random distributions of (a) LER 
{H = 5%) and (b) defects (n<j = 2.5%). 



to Ref. ifTsl . the effective electron mass ruj on the j-th 



subband reads: 



2 h^E, 



3 aH^A, 



(10) 



where t is the graphene hopping parameter (-2.7 eV) and Aj — 
cos (ttj/ (/ + 1)), where / is the number of dimer lines of the 
GNR. For the first conduction subband Ecjo — Eg/2, where 
Eg is the energy gap and j (which runs from 1 to I) is the 
index for which Aj is closest to -1/2. 

The corresponding Density of States (DOS), accounting for 
energy relaxation at outermost layers of the GNR [TSl, reads: 



[E + Ecjo-EcjY 



I Ecjo [E - Ecj ) {E+2Ecjo - Ecj ) 



(11) 



By means of Eqs. (|9|l and ( fTTT ). the phonon-limited mobility 
of a ID conductor (Eq. (|7]i) can be expressed as a sum over 
all contributing subbands j fTTJ: 



2q 



TThn2DWkBT 
' Ec jO 



r+co 

J2 dErp.iE) 



f{E) [l-f{E)] 



E 



{E^E, 



E, 



CjO) 



Ecjo 



Ecj + Ecjo 

1/2 



(12) 



where n2D = nm/W is the total 2D electron density, W the 
GNR width and T is the temperature. 

For what concerns longitudinal phonons, scattering rates 
are evaluated as in Ref. ||6]. According to Ref. [61, only 



intrasubband scattering has been considered. In particular, the 
longitudinal optical (LO) phonon scattering rate reads as 



l/ropiE) 



^piD3iE±huJLo) 



ApWujLO 
(1 + cos 6'k,k' 



l-f{E±n^Lo) 
1 - .f{E) 



, (13) 



where = l/[exp{hujLo/kBT) — 1] is the Bose-Einstein 
occupation factor and n+ = + 1, huLo is the optical 
phonon energy, Dqp is the optical deformation potential and 
p = 7.6 X 10^* g/cm^ is the 2D density of graphene. 
The factor (1 + cos ^k.k') arises from the spinor nature of 
the graphene eigenfunctions and 6'k,k' = 6j — 9ji, where 
dj = arctg {kx/kyj). Here kx {k'^) indicates the initial (final) 
longitudinal electron wavevector referred to the Dirac point, 
whereas kyj — 27rj/ [(/ + 1) a] and kyj' (which is equal to kyj 
for intrasubband scattering) are the quantized initial and final 
transverse wavevectors, respectively, where a is the graphene 
lattice constant, I is the number of dimer lines and j = 1, 
The intravalley longitudinal acoustic (LA) phonon scattering 
rate can be expressed as 



l/TACjiE)-- 



ApWvs 



PlD]{E) (1+cos ^k.k'J 



(14) 



where = n+ + n^, D^c is the deformation potential for 
acoustic phonons, t;^ = 2 x 10'' m/s is the sound velocity 
in graphene and \qx\ — 2\kx\ is the module of the phonon 
wavevector under the backscattering condition. 

For both acoustic and optical phonons, we have considered 
the four lowest subbands. The electron momentum relaxation 
time Tpj is computed by adding the relaxation rate due to 
electron scattering with acoustic and optical phonons [TTl . As 
a final remark, the effective mobility including all type of scat- 
tering sources has been extracted by means of Mathiessen's 
rule llntot = 'i^/f^LER + l/pd + '^/^J^iMP + '^/^J.ph, where 
Pd and Pimp are the defect and impurity limited mobilities, 
respectively. We have verified the validity of Mathiessen's 
rule considering samples with more sources of non-idealities 
(i.e. LER, ionized impurities and defects) at the same time. 
Then we have compared the computed mobility with that 
obtained by adding single contributions with Mathiessen's 
rule, observing a relative error smaller than 3%, which lies 
within the statistical error 



III. Results and Discussions 

The simulated segment is a double-gate GNR, embedded 
in Si02 with an oxide thickness tox of 2 nm, 10 nm-long 
(Fig. [3]i. The segment length has been chosen to satisfy 
the assumption of loss of phase coherence at the segment 
ends. Indeed, according to recent experiments fil8L the phase 
coherence length is close to 11 nm in graphene. From a 
computational point of view, different widths W have been 
considered, ranging from 1 to 10 nm: 1.12 nm, 2.62 nm, 4.86 
nm and 10.10 nm. All simulations have been performed at 
room temperature T — 300 K. 
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Fig. 3. 3D structure of the simulated GNR segment. 



A. Line-edge roughness limited mobility 

LER-limited mobility as a function of W for different edge- 
defect concentrations H is shown in Fig. in the above- 
threshold regime, for a 2D carrier density n2D of 9 x lO^'^ 
cm^^. As in all figures in the paper, the error bars represent 
the estimated root mean squared error of the average of 
the statistical sample (|5]l. 

As predicted by the analytical model in Ref. 0, fiLEB. 
scales as W^. Such behavior holds for large H (^20%) and 
narrow GNRs (W < 5 nm), when scattering from edge defects 
is expected to be heavier, while, for wider GNRs and for 
smaller H, such a law is not obeyed. In particular, for GNR 
width larger than 5 nm, fiLER tends to saturate, since the 
increasing number of subbands contributing to transport coun- 
terbalance the number of final states available for scattering, 
enhancing scattering rates. As shown in Fig. I^', in narrower 
GNRs, the higher the electron density, the larger the effective 
mobility, because of stronger screening, uler decreases for 
high n2D and wider GNRs, due to mode-mixing, as already 
observed in Silicon Nanowire FETs |fT9l . Indeed, for wider 
GNRs biased in the inversion regime, more transverse modes 
are able to propagate in the channel due to the reduced energy 
separation between different subbands. This leads edge defects 
to become a source of intermode scattering, thus reducing 

Fig. 13}; shows ^iler as a function of H, where ^ler. 
1/H for wide GNRs, consistent with the Drude model, and 
also observed in graphene in the presence of defects ||20) . 
However, as soon as W decreases, quantum localization 
becomes relevant fH], and the Anderson insulator-like be- 
havior in is recovered (pler oc in agreement with 
analytical predictions Q. 



B. Defect-limited mobility 

Defect-limited mobility is plotted in Fig. [5^ as a function 
of W for different defect concentrations. Even in this case, 
localization affects mobility in narrower ribbons, especially 
for higher (2.5%). 

For a fixed defect density, mobility slightly increases with 
electron density, due to the larger screening (Fig. |5j)) and. 
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Fig. 4. a) LER-limited mobility as a function of W for n2D = 0.9 X 
10^^ cm~^ and for different H. Data extracted from Ref. (6l are also reported, 
b) LER-limited mobility as a function of n2u for H = 5%. c) LER-limited 
mobility as a function of edge disorder concentration H for n2D = 0.9 X 
10^3 cm-2 and for different GNR width W. 
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Fig. 5. a) Defect-limited mobility as a function of W for n2D = 9 X 
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71213 for a defect fraction = 0.5%. c) Inverse of the mobility as a function 
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Carrier density n2D in the GNR channel for different GNR width: 
12 nm, 2.62 nm and 4.86 nm (from top to bottom). 
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Fig. 7. a) Impurity-limited mobility as a function of W for different impurity 
concentrations njMP and impurity charges, b) Mobility as a function of the 
impurity charge for nji^jp = 10^^ cm~^ and for different W . In a) and b) 
"2D = 0.9 X 10^^ cm~^, except otherwise specified, c) Impurity-limited 
mobility as a function of n2D for different W (nj^p = 10^^ cm~^ and 
the impurity charge is +0.4g). d) Inverse mobility as a function of ra/A/P 
for W = 10.10 nm. The carrier density is n2D = 0.9 X 10^^ cm~^. The 
experimental slope 2 X 10^^'' Vs extracted in Ref. 1241 is also reported. 



for larger GNRs biased in the inversion regime, it saturates 
with increasing W, for the same reason discussed above for 
LER scattering. In Fig. |5};, /i^ is plotted as a function of n^. 
The wider the ribbons the closer mobility follows the simple 
Drude model {p,d oc l/n^), as expected for strong disorder 
and uncorrelated scatterers in 2D graphene sheets [l22l . For 
W = 10.10 nm atomistic simulations are in agreement with 
experimental results: a linear curve fitting (fi = C/ud) leads to 
a proportionality factor of 2.23 x 10~^^ Vs, similar to those 
extracted in the case of Ne+ and He+ irradiated graphene 
samples (7.9x10-16 Vs and 9.3x10-^^ Vs, respectively) |20|. 

In Fig. |6l the GNR carrier density for widths ranging 
from 1.12 nm to 4.86 nm is shown. As can be seen, in 
Fig. |6] Anderson localization strongly degrades electron mo- 
bility |21|, creating percolating paths in wider GNRs and 
blocking conduction in the narrower ones. 



C. Ionized impurities limited mobility 

Impurity-limited mobility /iiMP, as a function of W, is 
shown in Fig. for n2D = 9 x 10^^ cm"^, and for 
different impurity charge concentrations. As can be noted, 
even a high impurity concentration of 10^^ cm"^ yields large 
mobility for QAq impurity charge. However, no indications 
are present in literature regarding the amount of unintentional 
doping charge |23|, [24|. Therefore, in order to check also 
the effect of impurity ionization on the electron transport, 
statistical simulations have been performed by increasing the 
impurity charge up to +2q. Mobility as a function of the 
impurity charge is plotted in Fig. |7j) for different W and 
for 77,213 = 9 x 10^^ cm^^. In this case smaller values of 
jjL (1700 cm^A's) are obtained for very narrow GNRs, due 
to the strongly nonlinear impact on screening in the channel. 
Even in this case localization strongly degrades mobility for 
narrower ribbons. 

To further test the importance of unintentional doping in 
limiting mobility, we have considered excess charge densities 
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Fig. 9. Mobility limited by phonons, LER, defect and impurity scattering 

in the inversion regime for a LER concentration H=5%, = 0.5%. The 

parameters for the scattering rates have been taken from [6 . The experimental 

mobility from [4J is also reported. 7120 = 0.9 X 10^'^ cm^^, tijmp = 10^^ 
_2 



up to cm^, which have been encountered in exper- 

iments li23l . As shown in Fig. |7^, in this case mobility 
decreases down to 10^ cm^/Vs for narrower GNRs. In Fig. |7]; 
impurity-limited mobility is plotted as a function of n2D for 
niMP = 10^^ cm"^ and impurity charge +0.4q and for 
different W. According to lf25l . Ii26l . /zimp in graphene does 
not depend on the electron density. The behavior is different 
in GNRs because up to an electron density of 10^^ cm"^, 
only the ground state is occupied, so that the Size Quantum 
Limit approximation is verified Ii27i . li6J. Since the scattering 
rate 1/t oc e^^ li25l and the static dielectric function e 
increases with n2D li27ll . the screening becomes stronger with 
increasing n2D- As a consequence, pimp oc t in GNRs has 
the increasing monotonic behavior shown in Fig.|7];. In Fig.|7}l, 
we compare experimental results available in literature ll24l for 
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graphene, showing the inverse of the impurity-limited mobility 
as a function of tijmp for W = 10.10 nm and by considering 
an impurity charge of +QAq: as expected for uncorrelated 
scatterers, /iiMP I/uimp and, as can be seen, experiments 
and simulations show quite good agreement. 

D. Acoustic and optical phonon-limited mobility 

Our study has also been directed towards the investigation of 
the impact of phonon scattering, through the Kubo-Greenwood 
formalism ["281, f29l. A wide range of phonon parameter 
values is currently present in the literature ||l], ll30l . |0, ||26l 
(i.e. acoustic {Dac) and optical {Dqp) deformation potentials, 
as well as optical phonon energy hw^o)- We observe that 
the most widely used phonon parameters are those adopted in 
Refs. ig, CI, E], i.e. Dac = 16 eV, Kujlo = 160 meV and 
Dqp = 1.4 X 10^ eV/cm, where TujOlo is the zone-boundary 
LO phonon energy. Such values have been tested towards those 
provided in Refs. Q, ll32l . Il33l . showing good agreement as 
far as mobility is concerned. 

To prove the validity of our approach, we have first com- 
pared our results with those obtained by means of an accurate 
2D Monte Carlo simulator [7J. For a fair comparison, the 
same phonon parameters and the same scattering rates as 
in Q have been used. As can be seen in Fig. |8^, results 
are in good agreement, especially for wider GNRs. How- 
ever, such parameters correspond to the out-of plane mode 
ZO which, according to symmetry -based considerations ll34l . 
density functional study [35 1 and Raman spectroscopy |f36l, is 
much weaker than in-plane vibrations. 

Therefore, in the following, we adopt the parameters dis- 
cussed above for LA and LO phonons and the scattering 
rates described in Eqs. ( fT3] l and (fl4t . In Fig. [S}) acoustic and 
optical phonon limited mobility is shown as a function of n2D- 
As expected, emission scattering rates are found to be larger 
than absorption scattering rates, due to their higher Bose- 
Einstein occupation numbers. In addition, as also observed in 
graphene [37], we have verified that the contribution of optical 
phonons is negligible also in GNRs, and jiph is dominated 
by (intravalley) acoustic phonon scattering [|6|, Q (Fig. [8j)). 
Note also that, unlike in grafene where \/n2D |l38l, 

in GNRs the transverse confinement leads to a non-monotonic 
n2_D -dependence as in CNTs ||39l . As can be seen, slightly 
increases due to the reduced number of available states for 
scattering. 

We observe that several recent studies [138] , ||40| , W\\ have 
demonstrated that surface phonons of the substrate represent 
a severe source of scattering, which strongly limits transport 
in graphene. However, we expect this effect to be much larger 
in high-fc dielectrics like Hf02, rather than in Si02, which 
is the insulator considered in this work. This issue will be 
the topic of a more comprehensive work on electron-phonon 
scattering in GNRs, which is beyond the scope of the present 
paper. 

Finally, we compare the total mobility with experiments 
from Wang et al [4] (Fig. |9]l. In particular, we show the 
mobility limited by different scattering mechanisms as well 



as the total mobility computed by means of Mathiessen's 
rule. As can be seen, when using parameters in Ref. ['51, 
LER is the most limiting mechanism {H — 5%) for very 
narrower GNRs, while for wider GNRs defect scattering is 
predominant, if a = 0.5% is considered. As an additional 
remark, we have checked that the same conclusion holds 
even if we consider much lower deformation potentials for 
phonons, that decrease the impact of phonon scattering, such 
as those those provided in Ref. I,30J . 

IV. Conclusion 

We have defined a simulation methodology based on atom- 
istic simulations on statistically significant ensembles of GNR 
segments to understand the functional dependence of GNR 
mobility upon different factors, and to quantitatively assess 
the importance of different scattering mechanisms. 

We used such methodology to investigate mobility in GNRs 
of width ranging from 1 to 10 nm. First, we find that, unlike 
in 2D graphene, electron-impurity scattering in GNRs is far 
too weak to affect low-field mobility. In addition, using well 
established parameters for electron-phonon coupling, we find 
that phonon scattering is hardly the limiting factor of GNR 
mobility. For narrower GNRs, line-edge roughness is the 
main scattering mechanism. This result is consistent with the 
findings in [jSj, where wider nanoribbons with very rough 
edges are characterized. Finally, for a fixed defect density 
or LER, mobility tends to decrease with the GNR width for 
narrower devices, suggesting the occurrence of localization 
effects. 
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